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ABSTRACT 



Larson (1981) first noted a scaling relation between masses and sizes in molecular clouds that implies that these objects 
have approximately constant column densities. This original claim, based upon millimeter observations of carbon 
monoxide lines, has been challenged by many theorists, arguing that the apparent constant column density observed is 
merely the result of the limited dynamic range of observations, and that in reality clouds have column density variations 
over two orders of magnitudes. In this letter we investigate a set of nearby molecular clouds with near-infrared excess 
methods, which guarantee very large dynamic ranges and robust column density measurements, to test the validity 
of Larson's third law. We verify that different clouds have almost identical average column densities above a given 
extinction threshold; this holds regardless of the extinction threshold, but the actual average surface mass density is a 
function of the specific threshold used. We show that a second version of Larson's third law, involving the mass-radius 
relation for single clouds and cores, does not hold in our sample, indicating that individual clouds are not objects that 
can be described by constant column density. Our results instead indicate that molecular clouds are characterized by a 
universal structure. Finally we point out that this universal structure can be linked to the log-normal nature of cloud 
column density distributions. 
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1. Introduction 

It has long been recognized that star formation is inex- 
tricably linked to the molecular clouds where the process 
is taking place, and therefore it is important to study the 
structure of these objects. One o f the first atte mpts in this 
direction has been carried out bv lLarso nl ((l98l . In his sem- 
inal work, Larson used molecular line data available from 
earlier studies (mostly millimeter data of nearly objects) 
and showed that molecular clouds obey three scaling rela- 
tions: (1) a power-law relationship between the length L of 
the cloud and its velocity dispersion er v , with a v oc L 38 ; 
(2) approximate virial equilibrium, with 2GM/a^L ~ 1; 
and (3) a relationship between the density n of the cloud 
and its length, with n oc L~ 11 . Larson's third law, which 
is the main focus of this letter, implies that molecular 
clouds have approximately constant column densities E, 
since £ ~ nL oc L~ 01 . 

Since their formulation, Larson's laws have been the 
subject of several observational an d theoretical s t udies . 
From the observational point of view, ISolomon et"aH (|1987|) 
presented 12 CO data for a 273 nearby clouds from the 
University of Massachusetts-Stony Brook (UMSB) Galactic 
Plane Survey (|Sanders et al.|[l986t ). They found a size- line 



with £(H2) = 170 M Q pc~ 2 . R ecently, the s ame s ample 
of clouds has been reanalysed by iHever et al.l (|2009f ) using 
data fro m the Boston Univ ersity-FCRAO Galactic Ring 
Survey (I Jackson et al.l 120061) . The use of 13 CO (J = 10) 
emission instead of 12 CO ensures that a large fraction of 
the data are optically thin; additionally, the data used have 
a much higher an gular sampling and spectral resolution. 
IHever et al.l (|2009t ) confirmed Larson's third law with a rel- 
ative scatter (approximately a factor 3) similar to previous 
studies. However, surprisingly they found a median mass 
surface density of molecular hydrogen for this sample of 
42 M(7) pc 2 , thus s i gnific antly smaller than the one derived 



bv ISolomon et al.l (jl987l ) 



On the theoretical side, there have been many at- 
tempts to explain Larson's laws using numerical simula- 
tions. In many cases, the validity of Larsons relations, 
and especial l y of t he third law, has been questioned 
(iKeeell Il989t IScalol [19901 IVazauez-Semadeni et al.l 119971: 
Ballesteros-Paredes fc Mac Low! 20021 : iBallesteros-ParedesI 
2006). In particular, it has been suggested that this law is 



width relationship with a power index (0.5) steeper than 
the one derived bv iLarsonl (|1981l ). Additionally, in agree- 
ment with Larson's third law, they found that the molecular 
gas surface density is approximately constant for all clouds 
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merely the result of the limited dynamic range of observa- 
tions, and that in reality mass surface densities of molecular 
clouds span at least two orders of magnitude. 

In this letter, we re-examine the validity of Larson's 
third law using e xtinction as a tracer of molecular gas 
(|Lada et al.l 119941) . The use of this tracer, in combina- 
tion with adva nced techniques (jLombardi fc Alvesl [2001; 
Lombardi 2009), allows us to probe clouds over a large dy- 
namical range (typically more than two order of magnitudes 
in extinction); additionally, the column density measure- 
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ments use a simple tracer, dust, which is not plagued by 
the uncertainties affecting millimeter observations of gas 
and dust (e.g., deviations from local thermodynamic equi- 
librium, opacity variations, chemical evolution, small-scale 
structure, depletion of molecules, unknown emissivity prop- 
erties of the dust, unknown dust temperature). 

The results of this study are twofold: first, we verify that 
Larson's law of constant column density holds with a very 
small sc atter on a set of nearby clouds investig ated using 
Nice r (jLombardi fc Alved l200lh and Nicest (jLombardil 
2009]); second, we show that the same law, applied within 
a single cloud (using different extinction thresholds) as 
M oc L 2 does not hold. Additionally, we argue that the 
first version of Larson's third law implies a universal physi- 
cal structure for molecular clouds, which we identify in their 
log-normal distributions for the projected gas density. 

Larson's third law, in its original formulation, links 
the average density ^n(H2)) of clouds with their size L: 
(n(H 2 )) = 3400 cm" 3 (L/T pc) Q , with a = -1.10. Here L 
is defined as the maximum observed linear extent of the 
cloud, and (n(H 2 )) is the average density of a sphere of 
diameter L and total mass M identical to the cloud (typi- 
cally estimated from 13 CO data). Larson's data were more 
heterogeneous and included different clouds studied at dif- 
ferent contours of integrated intensity, which resulted in 
a scatter of approximately one order of magnitude about 
the assumed relation; as we will see, our data suggest in- 
stead that Larson's law holds with a scatter below 15%. 
The fact that a ~ — 1 implies that the cloud projected 
column density, (n(H2))L oc L~ 01 , is approximately con- 
stant. Larson discussed a few possible explanations for this: 
one-dimensional shock compressions, optical depth natural 
selection effects, and observational biases owing to the lim- 
ited dynamic range of the 13 CO data. 

2. An extinction measurement of Larson's law 

2.1. Definitions 

We consider first (Sect. 12. 3p the following version of 
Larson's third law. Since we have at our disposal complete 
extinction maps, we can consider the area S of a cloud 
above a given extinction threshold Aq (throughout this let- 
ter, unless otherwise noted, we will refer to extinction mea- 
surements in the K band, Ak, and drop everywhere the 
index K). We then define the cloud size implicitly from 
S = 7r(L/2) 2 (or the cloud radius as R = L/2). Similarly, 
we can consider the cloud mass M above the same extinc- 
tion threshold. 

We will also briefly investigate the mass vs. radius re- 
lationship for each individual cloud, and verify whether we 
recover Larson's prediction M(R) oc R 2 (Sect. [23]). Note 
that the two versions of Larson's third law (different clouds 
above a fixed extinction threshold, or same cloud at var- 
ious extinction thresholds) are clearly linked, but are not 
equivalen t, in the sense that only one of the two might hold. 
Note also lLarsonl (|1981l ) de-facto studied different clouds at 
different thresholds, and therefore used a mixture of both 
versions considered separately here. 

Throughout this letter we will treat molecular com- 
plexes as single objects, and we will not split unconnected 
regions belonging to the same complex. Since typically a 
cloud will have many clumps with relatively high column 
densities, this procedure avoids the "creation" of new clouds 
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Fig. 1. Cloud masses above extinction thresholds of Ao = 
0.1 mag (filled symbols) and Aq = 0.5 mag (open sym- 
bols) as a function of their size. The two line shows the 
best constant surface density fits, which correspond to 
E = 41 M© pc~ 2 and E = 149 M pc" 2 respectively. 

when the extinction threshold Aq is increased. This proce- 
dure is justified because our objects are mainly well defined 
regions, relatively far from the galactic plane, and with no 
or little contamination from other clouds. 

2.2. Data analysis 

The data used here are extinction maps obtained from the 
point so urce catalog of the Tw o Micron All Sky Survey 
f2MASS: lKleinmann et~aTll 19941) . Data for the various com- 
plexes have been re duced using Nic er (|Lombardi fc Alvesl 
l200lh and Nicest (jLombardil I2009D and fo llowing the pre- 
scrip t ions adopted in previous works (see lLombardi et al.l 
120061 120081 I2010D . The complexes considered are nearby 
molecular clouds, and therefore we are able to well resolve 
most cores with the 2MASS d ata; the same clouds have 
been used in ILada et al.l (J2010). Extinction measurements 
are converted into surface mass densities using 

E = nm p l3 K A K , (1) 

where /j, is the mean molecular weight, (3k = [-^(Hi) + 
2N(K. 2 )]I Ak ~ 1.67 x 10 22 cm" 2 mag" 1 is the gas-to- 
dust ratio (jSavaee fc Mathislll979tlLillevlll955tlBohlin et all 
1978), and m p is the proton mass. With a standard gas 
composition (63% hydrogen, 36% helium, and 1% dust) we 
have fi ~ 1.37 and therefore E/Ar- ~ 183 M Q pc" 2 mag" 1 . 

2.3. Larson's third law for a constant extinction threshold 

Figure [T] shows the amount of mass different clouds have 
above extinction thresholds of Ak =0.1 mag and Ak = 
0.5 mag as a function of the cloud "radii" (defined accord- 
ing to Sect. 12. ip . together with the best power-law fit. As 
apparent from this plot, all clouds follow exquisitely well 
a Larson- type relationship, with M oc R 2 , and have there- 
fore very similar projected mass densities at each extinc- 
tion threshold. This result is also quantitatively shown in 
Table [TJ where we report the best-fit power-laws for the 
mass vs. radius relation at different extinction thresholds. 
The exceptionally small scatter observed in Fig. Q] is also 
confirmed by the results shown in Tabic [1] at all extinctions 
considered, data follow the best-fit power-laws with relative 
standard deviations always below 15%. 

Table Q] also show the dimensionless factor c obtained 
from the best quadratic fit M = c[im p l3KAoirR 2 in terms of 
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Fig. 2. Mass vs. radius relationship; both quantities are 
defined as indicated in Sect. 12.11 



the constants appearing in Eq. ((J}. Hence, c = (Ak)/Aq > 
1 , and the fact that c ~ 2 with a very small relative scatter 
among different clouds indicates that all these objects have 
a very similar physical structure. 

2.4. Larson's third law for single clouds 

Figure [5] shows the second version of Larson's third law con- 
sidered here, i.e. the mass vs. radius relationship. As appar- 
ent from this figure, the tracks for the various clouds have 
similar trends, but span a relatively large range of masses. 
In the range R € [0.1, 1] pc we can fit a power-law of the 
form M(R) = 380 M Q (R/pc) 1 6 a result that compare s 
well with the one obtained by iKauffmann et al.l ((2010), 
M(R) = 400 M Q (R/pc) 1 - 7 . Different clouds have quite 
similar exponents (the standard deviation of the power-law 
index is ~ 0.18), but rather different masses (the best-fit 
scale parameter for the mass ranges from 150 to 710 Mq). 
Note, however, that since the power-law index is signifi- 
cantly different from two, errors on the assumed distances 
of the clouds would affect the scale parameter for the mass. 

From this analysis we conclude that Larson's third law 
is not an accurate description of the mass vs. radius rela- 
tionship for single clouds. Specifically, at larger scales all 
clouds show a flattening of the curves and deviates signifi- 
cantly from a power-law, while at smaller scales clouds fol- 
low power-laws, but with an exponent significantly different 
than two. 
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1.63 
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1.44 


1.5 


379.8 


2.07 
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1.38 



Table 1. Best power-law fits M — airR 1 for various ex- 
tinction thresholds. Note that because 7 ~ 2 in all cases, 
the quantity a can be interpreted as the average mass col- 
umn density of the cloud above the corresponding extinc- 
tion threshold. The last two columns show the standard 
deviation of the cloud column densities divided by their 
average (relative scatter) and the ratio between the aver- 
age column densities and the minimum column density set 
by the extinction threshold (c). 
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Fig. 3. Cloud mass surface density above an extinction 
threshold as a function of the threshold, in logarithmic 
scale. The dotted line shows the relationship between the 
cloud column density in M pc~ 2 and the extinction in the 
K band. 



2.5. Cloud physical structure 

As mentioned earlier, that an ensemble of clouds satisfies 
Larson's third law at different extinction thresholds sug- 
gests that clouds have a universal physical structure. 

In order to investigate this point better, we consider in 
Fig. [3J the average column density of cloud material above 
a given extinction threshold, as a function of the extinc- 
tion threshold. Figure [3J indicates a remarkable uniformity 
among the various clouds: they all show a relatively flat 
plateau up to ~ 0.1 mag, and then a constant rise up to 
2-5 mag. In the range Aq s [0.1, 1] mag, the curves for all 
clouds are confined within a relatively narrow region. In this 
extinction range we can fit a simple power-law to the data 
plotted in Fig. [31 obtaining E = 265 M Q pc~ 2 (A /mag) - 8 . 
Note that an error analysis of the data points in Fig. |3J at 
Aq < 0.05 mag shows that they are significant, because 
the large number of independent measurements contribut- 
ing to these data make the statistical errors negligible, and 
because the flatness of the plateau at low extinction val- 
ues makes them robust with respect to systematic errors 
(such as offsets in the Nicer maps due to extinction in the 
control field). 



3. Theoretical interpretation 

The results presented above indicates that clouds have sim- 
ilar s tr uctures. Observati o nally (s ee, e.g., Lombardi et al.l 



2008 : iKainulainen et all 120091: iLombardi et al.l 12010: 



iFroebrich fc Rowlesll2010l ). many clouds show a log-normal 
distribution at low extinctions: 



Pa(A) = 



1 



2naA 



exp 



(InA-lnAi) 2 



2ct 2 



(2) 



where A\ and a are two positive parameters. A tail at 
high extinctions, present in many clouds, is generally as- 
sociated with the effects of gravitational instability. The 
log-normality of pa ( A) is often linked with supersonic tur- 
bulence, although recent results show that this is also a 
common feature o f very different classes of cloud models 
(jTassis et al.ll2010D . 

Interestingly, we can express the mass and the area of a 
cloud above an extinction threshold as simple integrals of 
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Fig. 4. Dimensionless column density k = 'E/Aifim p 0, de- 
fined in Eq. ([5]), as a function of the dimensionless column 
density a = A /Ai. The dotted line shows the asymptotic 
limit «(a) ~ a. Note the similarity with the curves shown 
in Fig. [31 



Pa(A). Given a cloud of total area Stot, the area and mass 
above a given extinction threshold Aq are 



S(A ) = S tot Pa(A) dA 



M(A ) = Stot/impP / A PA (A) dA . 
Ja 



(3) 
(4) 



In particular, if we consider the log-normal distribution of 
Eq. ([2]), we obtain for the column density above A 



where K is a dimensionless quantity defined as 



k(o) 



exply 



cr 2 \ 1 - erf [(In a -<j 2 )/V2<j] 



1 - erf[lna/V2a] 



(5) 



(6) 



We plot in Fig. [4] the function n(a) in a log-log scale for 
three values of the parameter a. A comparison of Fig.Uwith 
Fig. [3J shows that the log-normal model is able to capture 
the main characteristics, except the "collapse" of curves 
onto the dotted line in Fig. |3J (most likely due to limited 
dynamic range of observations). Furthermore, in order to 
obtain the narrow bundle of curves in Fig. [3j the parame- 
ters A\ and a for the various clouds must span a limited 
range. Typical relative scatters of A\ and a are of the order 
of 57% and 19% (cf. Tab. [2]). Additionally, an analytical cal- 
culation shows that the particular form of the log-normal 
distribution further suppresses these scatters (respectively 
by factors between 8 and 4), so that the final expected 
relative standard deviation in E of the order of 14%, in 
agreement with the data presented in Table [TJ 

4. Summary 

1. Using near- infrared extinction maps of a set of nearby 
clouds we tested Larson's third law for molecular clouds, 
the constancy of average mass surface densities above 
a given extinction threshold. We verified this scaling 
law to a relatively high degree of precision. We found 
a very small (< 15%) relative scatter for the measured 
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Coalsack 
Taurus 
Lupus I 
Ophiuchus 
Serpens 
Cha I 
Cha II 
Lupus III 
Cor A 
LDN 1228 
Pipe 



0.15 
0.38 
0.18 
0.08 
0.16 
0.33 
0.11 
0.12 
0.14 
0.10 
0.12 
0.42 



0.42 
0.28 
0.49 
0.43 
0.48 
0.51 
0.35 
0.35 
0.35 
0.44 
0.32 
0.29 



LDN 134 
LDN 204 
LDN 1333 
LDN 1719 
Musca 
Cha 3 

Ori A CMC 
Perseus 
Ori B GMC 
Cepheus 
California 



0.10 
0.14 
0.12 
0.15 
0.08 
0.12 
0.13 
0.13 
0.11 
0.16 
0.14 



0.39 
0.41 
0.38 
0.50 
0.45 
0.46 
0.50 
0.48 
0.49 
0.59 
0.51 



Table 2. Lo g- normal fit parameters for various clouds. 
Data are from lKainulainen et al.l (|2009l ) and converted into 
if -b and extinction paramete rs using a standard reddening 
law dRieke fc Lebofskvl[T985h . 



column densities independent of the adopted extinction 
thresholds over a very large range, from Ak = 0.1 mag 
to Ak — 1.5 mag. Additionally, we found the value of 
the average mass surface density to be a function of the 
adopted extinction threshold. 

2. We verified that Larson's third law does not hold 
when considering the mass-radius relation within sin- 
gle clouds. In the range R e [0.1, 1] pc we find that the 
mass scales as M(R) cx R , and is therefore signifi- 
cantly shallower than what was predicted by Larson; at 
larger radii, the relation appears to flatten even more. 

3. We interpreted these results, and in particular item 1 
above, as the effects of a universal physical structure 
shared among the different clouds. We showed that this 
universal structure is represented by a uniformity in the 
cloud density distributions. We found that a log-normal 
model is able to account for this uniformity, provided 
that the log-normal parameters are restricted to rela- 
tively narrow ranges. This suggests that Larson's third 
law might be a consequence of this special property of 
cloud structure. 

Acknowledgements. We thank the referee for helping us to signifi- 
cantly improve this paper. 
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